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A four-component Fermi gas in one dimension with a short-range four-body interaction is shown 
to exhibit a one-dimensional analog of the BCS-BEC crossover. Its low-energy physics is governed 
by a Tomonaga-Luttinger liquid with three spin gaps. The spin gaps are exponentially small in the 
weak coupling (BCS) limit when it is due to the pairing, while they are large in the strong coupling 
(BEC) limit because of the formation of tightly-bound tetramers. We investigate the sound velocity 
and the gap spectrum in the BCS-BEC crossover and discuss exact relationships valid in our system. 
We also show that a one-dimensional analog of the Efimov effect occurs for five bosons while it is 
absent for fermions. 
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I. INTRODUCTION 



II. FEW-BODY PROBLEMS 



A. Lattice model 
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Experiments using ultracold atomic gases have 
achieved striking progress in realizing and studying var- 
ious many-body systems previously regarded as theoret- 
ical models. One example is the two-component Fermi 
gas with a short-range two-body interaction that exhibits 
the BCS-BEC crossover [1, 2] . In the weak coupling limit, 
the system is a BCS superfluid where fermionic excita- 
tions have an exponentially small gap, while at strong 
coupling it becomes a dilute Bose-Einstein condensate of 
tightly-bound dimers with a large gap for the fermionic 
excitations. These two limits are smoothly connected by 
varying a single parameter, the scattering length. When 
the scattering length is much larger than the range of 
the interaction potential, the properties of such a sys- 
tem become independent of the shape of the interaction 
potential. This universality makes the study of the BCS- 
BEC crossover extremely worthwhile because the same 
properties can be shared by many different systems [3]. 



In this paper, we investigate a purely one-dimensional 
analog of the BCS-BEC crossover in a four-component 
Fermi gas with a short-range four-body interaction. The 
short-range four-body interaction in one dimension is 
characterized by the scattering length exactly in the same 
way that it characterizes the short-range two-body inter- 
action in three dimensions [4] , and therefore, leads to the 
universal "BCS-BEC" crossover in one dimension. We 
note that the BCS-BEC crossover of a two-component 
Fermi gas in a quasi-one-dimensional geometry has been 
studied in Refs. [5, 6], while the crossover studied in this 
paper is universal, being independent of the confinement 
potential. We also note that four-component (spin-3/2) 
Fermi gases with two-body interactions have been stud- 
ied and reviewed in Ref . [71 . 



We start with a system of fermions with four compo- 
nents labeled by cr = a, b, c, d living on a one-dimensional 
lattice. We assume that each lattice site can accommo- 
date one, two, or three particles with no change in energy, 
but an introduction of a fourth particle into a site with 
three particles gains a finite amount of energy. The lat- 
tice Hamiltonian for such a system is given by 



H = -t 

(xy),a X 

(1) 

We will be interested in the dilute limit where the average 
number of particles per site is small. In order to find 
the universal regime, we consider the scattering of four 
particles with all different components. Such a four-body 
problem is described by the Schrodinger equation 



*(a;) =E'^{x), 



(2) 



where x — {xa,Xb,Xc,Xd) is a set of coordinates of four 
particles and is the discrete Laplacian with respect 
to x„; A„^{x„) = ^{x„ + 1) + ^(x^ - 1) - 2*(a;^). The 
four-body interaction potential is given by V{x) — —U 
when all equal and V{x) — otherwise. 

Since V{x) is translationally invariant, it is convenient 
to introduce new coordinates X — {xa + Xh + Xc + Xd)/4:, 

ri = (Xa + Xb - Xc - Xd)/2, r2 = {Xa - Xh + Xc - Xd)/2, 

and ra — (xa — Xb ~ Xc + Xd)/2 and assume ^'(a;) to 
be independent of the center-of-mass coordinate X. The 
Schrodinger equation (2) in terms of the remaining three 
relative coordinates r — (ri,r2,7'3) becomes 



*(r) 



where Ai4'(r) 
^(1,1,1) 



ei 



62 



^{r - e 
-1,1), 63 



£;*(r) 
) 



(3) 



2^'(r) with 
i(-l,l,-l). 
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and 64 = 1,— 1). Eq. (3) is equivalent to the 

Schrodinger equation describing one particle moving in 
a body-centered cubic lattice with an attractive poten- 
tial of the magnitude U concentrated at one lattice point. 

One can see from Eq. (3) that the zero-energy wave 
function at a long distance has the form 



\l/(|r| oo)\ 
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(4) 



which is familiar in a two-body scattering in three di- 
mensions. This can be understood from the fact that 
the continuum limit of Eq. (3) with t = l/(2m) exactly 
resembles the Schrodinger equation in three dimensions. 
Here a is an arbitrary real parameter characterizing the 
long-distance physics and referred to as the scattering 
length. By matching the solution of Eq. (3): 
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with the asymptotic form (4), we find a in units of the 
lattice spacing to be 
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The scattering length a can be fine-tuned to infinite cor- 
responding to the four-body resonance by choosing 
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5.742. 



(7) 



This value oi U/t separates the weak coupling regime 
(a < 0) with no bound state from the strong coupling 
regime (a > 0) with a four-body bound state (tetramer) 
whose wave function and binding energy for a ^ 1 are 
given by the universal formulas: 



*(|r| 00) cx 
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The long-distance physics near the critical value of U/t 
should be universal and, in particular, the scale and 
conformal invariance are achieved in the unitarity limit 
a —i' 00 [41 . 



B. Field-theoretical formulation 






FIG. 1: Feynman diagrams describing the four-body scatter- 
ing in vacuum. The dot represents the bare vertex ico and 
the dashed line represents the scattering amplitude iA{E,p). 



Throughout this paper, we neglect two-body and three- 
body interactions and interactions involving the same 
components of fermions. In addition to the translational 
and Galilean symmetries, the Hamiltonian density has 
global U(l) and SU(4) symmetries under 



and 



(10) 



corresponding to the conservations of charge and spins, 
respectively. 

The second term in Eq. (9) describes the four-body 
contact interaction among all different components of the 
fermionic field ipa ■ co is a cutoff-dependent coupling con- 
stant and can be related to the above-introduced scatter- 
ing length a by matching the property in the four-body 
problem. The four-body scattering amplitude A{E,p) is 
obtained by the geometric summation of Feynman dia- 
grams (Fig. 1): 

.... .,1 1 . • /■°° dkidk2dk3 

2m 

^ kl + k'^ + k1 + (fci +k2 + ksf -f pV4 - 2mE - iO+ ' 

(11) 

Here the integrations over momenta fci, fc2, and fca 
are linearly divergent. Introducing a momentum cutoff 
y^kf + k2 + fc| < A and choosing the cutoff dependence 
of Co as 



Co 



mA 



(12) 



we obtain the following cutoff-independent scattering 
amplitude in the limit A — >■ cx3: 



AiE,p) 



4:77 
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-1/a + 0?74^^2^^IE^^^ 



(13) 



In particular, when a > 0, A{E,0) has a pole on the 
negative real axis of E indicating the existence of the 
four-body bound state. Because its binding energy is 
given by Eq = —l/i/Zma^), we can identify a in Eq. (12) 
with the scattering length introduced in Eq. (4). 



The physics in the universal regime can be described 
by the following Hamiltonian density in the continuum 
limit: 



n 



^ 2m 



Co i'l tpl tpltpdIpcMa ■ (9) 



C. Five-body problem and Efimov effect 

The above arguments equally apply to four-component 
bosons in one dimension. However, the many-body sys- 
tem of bosons with the attraction interaction tends to 
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For identical bosons with the four-body resonant interac- 
tion, the anomalous dimension of the five-body composite 
operator satisfies Eq. (15) with A = 4. It has com- 
plex solutions 7 — ±1.25 i, and therefore, the spectrum 
of pentamers is much denser; 



FIG. 2: Feynman diagrams to renormalize five-body compos- 
ite operators. The dashed line is a resummed propagator of (p 
field which is equal to —iA. The shaded bulb represents the 
vertex function z{p). 



be unstable, while the stability of fermions is feasible be- 
cause of the Pauli exclusion principle. Such a difference 
can be seen already in a five-body problem: Five bosons 
develop deeply bound states while fermions do not. 

In order to see this, we study a scaling dimension of 
five-body composite operator (jjip^ in the unitarity limit 
a — cx), where (j) = CQipaipb^c^d is a tetramer field. Feyn- 
man diagrams to renormalize (jjip^ is depicted in Fig. 2. 
The vertex function z{p) satisfies the following integral 
equation: 



z{p) + X 



J-oo2vr 



dkidk2 



(14) 



where A = ±1 for four-component bosons/fermions. Be- 
cause of the scale invariance, we can assume z(p) cx 
(|p|/A)^, where the anomalous dimension 7 satisfies 



4;^ cos(7arctan--^^ 



(15) 



The anomalous dimensions of even-parity operators sat- 
isfy the same equation (15) and their scaling dimensions 
are given by 



■7=2+7. 



(16) 



For fermions (A = —1), we can find a series of real 

solutions; 7 = 1.59,4.08,5.99, According to the 

operator-state correspondence [8-10], each solution cor- 
responds to the energy of resonantly-interacting five 
fermions in a one-dimensional harmonic potential by 



E : 



2+7)'^. 



(17) 



On the other hand, for bosons (A = -fl), in addition to 
real solutions 7 — 2.25, 3.91, 6.01, . . . , we can find a pair 
of complex solutions 7 = ±0.735 i. This is a signal of 
the formation of an infinite set of five-body bound states 
(pentamers) whose spectrum exhibits the discrete scaling 
symmetry: 
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This is an analog of the Efimov effect in three dimen- 
sions for three identical bosons [11]. For comparison, the 
scaling factor for Efimov trimers in three dimensions is 
known to be e^'^/l^™''! = (22. 7)^. 

Similarly, the anomalous dimensions of odd-parity op- 
erators [e.g. 4(/)(VV'o-) — (V0)?/'<t] satisfy 



4 A sin(7arctan^) 
7T5 '^"^ 



7Cos(^) 



(20) 



with the scaling dimensions given by Eq. (16). In this 
channel, both fermions and bosons have real solutions 
only; 7 = 0.833,3.15,4.87,... for A = -1 and 7 = 
1.17, 2.85, 5.12,... for A = +1. Therefore, the corre- 
sponding states in a harmonic potential are universal. 



III. MANY-BODY PROBLEMS 

Since bosons with the four-body resonant interaction 
in one dimension develop deep five-body bound states, 
the corresponding many-body system cannot be stable 
toward collapse. Therefore, we will study the many-body 
physics of four-component Fermi gas in one dimension 
as a function of the dimensionless parameter kpa char- 
acterizing the short-range four-body interaction. Here 
hp = Tin /A is the Fermi momentum defined by the to- 
tal number density n. According to the corresponding 
BCS-BEC crossover in three dimensions, it is natural to 
refer to the weak (strong) coupling limit kpa — (+)0 
as the "BCS" ("BEG") limit ahhough we do not have 
the spontaneous symmetry breaking in one dimension. 



A. BCS limit 

The many-body physics is conveniently described by 
introducing the chemical potential term —fxip^ipa to the 
Hamiltonian density in Eq. (9). In the weak coupling 
(BCS) limit fcpa — > —0, the system develops two Fermi 
points at p = ±fcF and low-energy degrees of freedom 
are excitations around the Fermi points. Expanding the 
fermionic field as 



1pa{x) 



i;f{x)+e-'''-^ij^{x), (21) 
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the low-energy effective theory consistent with the origi- 
nal symmetries (10) can be written as 

^BCS = -iVF V'?^ VV-^ + IVF VV-^ 

+ I (^^t^Rt^R^R + V^t^Lt^L^L^ ^ 

where = fcp /"^ is the Fermi velocity and summa- 
tions over (7(t) = a, h, c, d are implicitly understood. The 
low-energy parameters g\, g^, and 54 are determined by 
matching two-body scattering amplitudes at the Fermi 
points with those from the microscopic theory (9). To 
the leading order in kpa, we find 



9i = 92 



■ 94, = KF|a| 



0[{kFaf 



(23) 



The spectrum of the low-energy effective theory "Hbcs 
can be obtained exactly in terms of the bosonization. We 
introduce charge current operators 



J, 



R(L) _ 







^R(L)t^R(L) 



and spin current operators 
J^(L)^^R(L)t 



;R(L) 



(24) 



(25) 



where ta with a = 1, . . . , 15 are generators of SU(4) Lie 
algebra normalized as Tr(tQ,i/3) =<5a/3/2. Using these cur- 
rent operators, Hbcs <2an be separated into two mutually 
commuting parts [12]; 'Hecs = '^ch + 'Hgp with 



^ch = 



2'KV-p -\- 3(74 



tR tR I jLt JL 

Jq Jq -r Jq Jq 



L tL\ 



4.92 - 91 jR jL 



and 



(26) 
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2ttvf - 54 



' tR tR I 7L tL\ _ 9 tR tL 



(27) 

The charge part "Hch is easily diagonalized by 
the Bogoliubov transformation and equivalent to the 
Tomonaga-Luttinger liquid. The standard bosonization 
rule leads to the Hamiltonian density describing a gapless 
excitation transporting a particle number [12]: 



ch 



2ttK 



(28) 



Here the Tomonaga-Luttinger parameter K and the 
sound velocity Vg in the BCS limit kpa —0 are given 
by 



6A;f|ci| 



and 



6A:F|a| 



Vp. 



(29) 



We note that the relationship Kvg = vp is guaranteed by 
the Galilean invariance [13]. 

On the other hand, the coupling gi < in the spin part 
■Hsp is marginally relevant and thus develops gaps in the 



spectrum. The exact gap spectrum can be obtained from 
the Bethe-ansatz solution if we recognize T-Lsp as the non- 
Abelian part of the SU(4) chiral Gross-Neveu model [14]: 



Afocvphpe " '^-"'i-i'sm 



(30) 



Here / = 1, 2, 3 is a fermion number of excitations, and 

accordingly, there are three distinct gaps which arc ex- 
ponentially small in the BCS limit k-pa —0. The 
degeneracy of Ai and A3 can be traced back to an 
accidental symmetry in "Hecs under the charge con- 
jugation; ^CT^^'' ^ ip^^^'^^ . This symmetry is broken 
by quadratic derivative terms '^ct^^'^V^Vct^^VC^"^) ne- 
glected in Hbcs- Therefore, the degeneracy has to be 
slightly lifted by 



A3 — Ai vpkpe 



V(4fcF|o|) 



(31) 



BEC limit 



In the strong coupling (BEC) limit kpa — >■ -|-0, four 
fcrmions with all different components form a tightly- 
bound tetramer and thus the many-body system will be a 
dilute Bose gas of such tetramers. In this limit, fermionic 
excitations are largely gapped because of the binding en- 
ergy of the tetramer Eq = — l/(2ma^). The gap spec- 
trum with the fermion number / = 1, 2, 3 is simply given 

by 



/ 



Sma? 



(32) 



Interestingly, we find that the ordering of the three dis- 
tinct gaps is Ai < A2 < A3 in the BEC limit while it is 
Ai ~ A3 < A2 in the BCS limit [see Eqs. (30) and (31)]. 
Therefore, there has to be a crossing between two gaps 
A2 and A3 as a function of —00 < {kpa)~^ < 00 in the 
BCS-BEC crossover. 

The dilute Bose gas of tetramers is described by the 
Hamiltonian density 



1 



BEC 



2M 



Matt 



(33) 



where M = 4m is the tetramer mass and the tetramer 
density is n/4. att is a tetramer-tetramer scattering 
length (analogous to the dimer-dimer scattering length 
in three dimensions [15]) characterizing the scattering of 
two tetramers in one dimension. Because the scattering 
length a is the only scale of the system in vacuum, att 
has to be proportional to a; 



att = -m- 



(34) 



The coefficient 77 is a universal number obtained by solv- 
ing the eight-body problem of fcrmions nonperturba- 
tively. 77 is expected to be positive because tetramers 
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should repel each other due to the fermionic statistics of 
the constituents. If ry > 0, then the many-body system 
of tetramers is stable. Here we shall assume 77 > and 
leave the determination of the exact value of 77 as a future 
problem. 

The effective theory of tetramers "Hbec is nothing but 
bosons with a (5-function interaction in one dimension and 
can be solved exactly [16]. In particular, its low-energy 
physics is described by the Tomonaga-Luttinger liquid 
(28) again [13]. The Tomonaga-Luttinger parameter K 
and the sound velocity in the BEC limit k^a — >■ -|-0 
are given by 



) 



and I 1 — 



- J 4' 

(35) 

In the expression for K, we have taken into account the 
fact that the particle number of a boson is four. The 
BCS-BEC crossover indicates that K and Vg in Eqs. (29) 
and (35) are smoothly connected as functions of —00 < 
(fcpa)"^ < 00. 



C. Unitarity limit 

It would be difficult to compute K and Vg away from 
the BCS or BEC limit. However, in the unitarity limit 
/cpci — >■ 00, we can derive exact relationships between K 
and thermodynamic quantities. Because the density n 
is the only scale of the system, the ground state energy 
density of the unitary Fermi gas can be written as 



(36) 



where the ground state energy density of a noninteracting 
Fermi gas is 



96m 



(37) 



Here ^ is a universal number to characterize the strongly- 
interacting unitary Fermi gas and analogous to the 
Bcrtsch parameter in three dimensions [3]. From the 
thermodynamic relationships, we obtain the pressure as 
P(n) = 25(n), and thus, the sound velocity is given by 



1 dP 
m dn 



(38) 



Because Kvs = vp is guaranteed by the Galilean invari- 
ance [13], we find that the Tomonaga-Luttinger parame- 
ter is related to the one-dimensional Bertsch parameter 
by 



1 

VI' 



(39) 



This relationship implies if > 1 in the unitarity limit be- 
cause < ^ < 1 is expected for the attractive interaction. 
It is a challenging many-body problem to determine the 
exact value of 



D. Exact relationships 

Another characteristic of our system (9) that resem- 
bles the BCS-BEC crossover in three dimensions is the 
large-momentum tail of the momentum distribution of 
fermions and its relationships to other properties of the 
system [17]. Here we derive the exact relationship be- 
tween the coefficient of the large-momentum tail C and 
the energy of the system E. We consider the following 
operator product expansion (no sum over a = a,b, c, d): 



VM 

87r 



(40) 



This can be confirmed by evaluating expectation values of 
the both sides for a state consisting of four fermions with 
all different components [18]. The nonanalytic term ~ |y| 
indicates that the momentum distribution of fermions 



o,{k) = jdxj dye-'^y (v4 (.t - |) Vv [x + | 



falls off by a power of |fc| —5 
Pa{k) 

The coefficient is given by 



00 as 



C 



(42) 



C = 



47r 



dx {mcof{i)lipliplijlipdi^ci^btpa){x) (43) 



and called an integrated contact density. From Eqs. (9), 
(12), and (43), wc find that the energy of the system 
E = Jdx {H) can be expressed by 



Paik) 



,27r 2m 



fc2 



C 



(44) 



This relationship is valid for any state of the system and 
for any value of the scattering length a. Derivations of 
the pressure, the adiabatic relationship, and the virial 
theorem [17, 18] are straightforward from above results. 



IV. CONCLUSIONS 

In summary, we have demonstrated that the four- 
component Fermi gas in one dimension exhibits the one- 
dimensional analog of the BCS-BEC crossover as a func- 
tion of the scattering length characterizing the short- 
range four-body interaction. We investigated the sound 
velocity, the gap spectrum, and the exact relationships in 
the BCS-BEC crossover and foimd that the gap spectrum 
has the rich structure because of the existence of three 
distinct gaps. We also showed that the one-dimensional 
analog of the Efimov effect occurs for five bosons while 
it is absent for fermions. 
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